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The mathematical properties of the new analytic running coupling (NARC) in QCD 
are investigated. This running coupling naturally arises under "analytization" of the 
renormalization group equation. One of the crucial points in our consideration is the 
relation established between the NARC and its inverse function. The latter is expressed 
in terms of the so-called Lambert W function. This relation enables one to present 
explicitly the NARC in the renorminvariant form and to derive the corresponding j3 
function. The asymptotic behavior of this /3 function is examined. The consistent 
estimation of the parameter Aqcd is given. 

1. Introduction 

The description of hadron interactions at small characteristic momenta remains 
an actual problem of elementary particle theory. The asymptotic freedom in quan- 
tum chromodynamics (QCD) enables one to investigate processes at large momen- 
tum transfers by making use of standard perturbation theory. But there is a num- 
ber of phenomena which description lies beyond such calculations, namely quark 
confinement, gluon and quark condensates and many others. For these purposes 
nonperturbative approaches are used. 

In the late 50 's the so-called analytic approach to quantum electrodynamics 
was proposed.Eli3 Its basic idea is the explicit imposition of the causality condition 
which implies the requirement of the analyticity in the variable for the relevant 
physical quantities. Recently this approach has been extended to QCD.i This leads 
to the following essential advantages: absence of unphysical singularities at any loop 
level, stability in the infrared (IR) region, stability with respect to loop corrections, 
and extremely weak scheme dependence. The analytic approach has been applied 
successfully to sucIUQCD problems as the r lepton decays, e+e^-annihilation into 
hadrons, sum rules. □ 

Recently thcLajialytic approach has been employed to the renormalization group 
(RG) equation.ErQ The analyticity requirement was imposed on the RG equation 
itself, before deriving its solution. Solving the RG equation, "analytized" in the 
above-mentioned way (i.e., the proper analytical properties of the RG equation £is 
a whole have been recovered), one arrives at the new analytic running couplingp 
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(NARC). An essential point, that plays a crucial role in description of a number 
of nonperturbative phenomena in our approach, is the IR enhancement of the new 
analytic running coupling at ~ 0. It should be stressed here that such a behavz 
ior of the invariant charge is in agreement with the Schwinger-Dyson equations Jj 
and, as it was demonstrated recently,EH3 provides description of quark confinement 
without invoking any additional assumptions. 

The objective of this paper is to study in details the mathematical properties of 
the new analytic running coupling in QCD.i This enables one to present the NARC 
in the explicitly renorminvariant form and to construct manifestly the correspond- 
ing /? function. We restrict ourself to the one-loop level calculations, taking into 
account the following consideration. Due to the higher loop stability!! all the essen- 
tial features of the new analytic running coupling become apparent at the one-loop 
level already. The account of higher-loop corrections does not lead to any qualita- 
tive changes. At the same time, at the one-loop level one able to carry out all the 
calculations in an explicit form, that is one of the self-evident merits of the current 
consideration. 

The layout of the paper is as follows. In Sec. 2 the new analytic running cou- 
pling is briefly discussed. In Sec. 3 the function N(a), that plays a key role in our 
analysis, is introduced and its properties are investigated. In turn, the function 
N{a) is expressed in terms of the so-called Lambert W function. Proceeding from 
this, in Sec. 4 the one-loop new analytic running coupling is presented manifestly 
in the renorminvariant form and its basic properties are revealed. Further, the 
corresponding (3 function is derived in an explicit form and its properties are exam- 
ined. The consistent estimation of the parameter Aqcd is given. In the Conclusion 
(Sec. 5) the obtained results are formulated in a compact way, and further studies 
in this approach are outlined. 

2. New analytic running coupling in QCD 

In the analytic approach, extended to QCD by Shirkov and Solovtsov,i the 
basic idea is the explicit imposing of the causality condition, which implies the 
requirement of the analyticity in the variable for the relevant physical quantities. 
Later this idea was applied to the "analytization" of the perturbative series (i.e., 
the recovering of the nroper analyticity properties of this series) when calculating 
the QCD observables.Q 

As known, the Q^-evolution of QCD observables is usually described in the 
framework of the RG approach. However, the standard RG calculations involve, 
as a rule, the use of the perturbative expansions. As a result, the RG equation at 
any given loop level does not meet the analyticity requirement.il! In order to avoid 
such situation, the procedure of analytization of the RG equation as a whole has 
been elaborated.ifl This procedure implies the recovering of the correct analytic 
properties of the RG equation itself, before deriving its solution. As a result, the 
solution of the analytized RG equation is expressed in terms of a new analytic 



running coupling. At the one-loop level this running coupling reads 
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where Po — 11 — 2 nf/3 is the first coefScient of the P function, nt is the number of 
active quarks. At the higher loop levels only the integral representation for ^a^niQ^) 
has been obtained.tl So, at the k-loop level we have 
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where 7^W(c^) = i^Triy^T, 
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is the spectral density, as is the /c-loop perturbative running coupling, a 
a/3o/(47r), and zq = Q\l !s? is the normalization point. 
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Fig. 1. The normalized new analytic running coupling a{z) = '^OanCQ^) /'*cran(Qo) one-, 
two-, and three-loop levels. The normalization point is Qq = lOA^, z = Q'^ / K"^ . 

It is worth making here a short historical comment. When considering the ana- 
lytic aporoach to quantum electrodynamics, it was noted that the method proposed 
in Refs.B@ is not unique. For example, in paper§ the ambiguity of that method was 
shown explicitly. Furthermore, it was stressed that for avoiding the ambiguity of 
the approach developed in Refs.Bi one should impose an additional condition,0 in 
particular, by making use of the equationsof motion. Specifically, following this 
way the boson propagator has been derivedllj in which the singularity was removed 



in a multiplicative way (i.e., similar to EkL (|^)). Moreover, another common feature 
of our consideration and that of papertj is the following. Unlike the method,Hi3 
in our approach one is also able to catch the nonperturbative contributions to the 
spectral density for the running coupling (see also Eq. (p3|)). 

Figure |l| presents the new analytic running coupling calculated at the one-, two- 
and three-loop levels. These curves show that the analytic running coupling (|l|) 
possesses the higher loop stability. Moreover, it can be demonstrated that the 
singularity of the NARC at the point = is of the universal type at any loop 
level.! 



3. The Lambert W function 

For the representation of the new analytic running coupling (^ in the renormin- 
variant form and for the construction of the relevant /3 function, it proves to be 
convenient to use the so-called Lambert W function. As long ago as the miidle 
of 18th century this function is being employed in diverse physical problems.tll In 
current researches in O CD the interest to this function arose just a few years ago. 
So, it was revealedll3'tj that the exact solution to the perturbative RG equation for 
the invariant charge at the two-loop approximation can be expressed in terms of 
the Lambert W function. At the three-loop level, when the Fade approximation for 
the P function is used, the sohition to this equation can be expressed in terms of 
the Lambert W function also.EJ But in this case one has to take into account that 
the application of the Fade approximation to the /3 function drastically changes the 
behavior of perturbative running coupling at small Q^. 

The Lambert W function is defined as a many-valued function W{x), that 
satisfies the equation 

W{x) exp [Wix)] = X. (3) 



Only two real branches of this function, the principal branch Wq (x) and the branch 
W-^i{x) (see Fig. will be used in our consideration. The other branches of the 
Lambert W function take an imaginary values. One can show that for the branches 
Wq and VF-i the following expansions hold: 

Wa{e) = e-e^ + 0{e^), e ^ 0, (4) 
W-i{~e) = lne + 0(ln|ln£|), e-^0+, (5) 

Wo(^-^+e^ = -l + V2^ + 0{e), e ^ 0+, (6) 

W-i(^-^+e^ = -l-V2^+0{e), e ^ 0+. (7) 

Details concerning the mathematical properties of the Lambert W function can be 
found in the review Jld 
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Fig. 2. The function xe^ (dot-dashed curve) and two real branches, Wo{x) (solid curve) and 
W--i{x) (dashed curve) of its inverse function. 



In order to represent the NARC in a renorminvariant form and to derive the 
P function corresponding to ^a^^,^ one has to solve the equation (see Sec. 4) 
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with respect to the variable z. Let us multiply this equation through by the factor 
Inz, 

lnz=--6, z > 0, 6=--<0. (9) 
z a 

The equation obtained has a trivial solution z — 1 that does not satisfy initial Eq. (^) 

when a 7^ 1. Obviously, this is a consequence of the multiplication of Eq. (H) by 

Inz. Therefore, when solving Eq. (^) for the variable z one has to discard its trivial 

solution z = 1. Next, we can represent Eq. (0) in the form 



— exp I - ) ^ be 



(10) 



Taking into account the definition (^), solution to Eq. ( p^ ) can be expressed in 
terms of the Lambert W function 
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The branch index k of the W function will be specified below. 

In the physically relevant range a > the argument of the Lambert W function 
in Eq. (O) takes the values — 1/e < be'' < 0. The only two real branches, the 



principle branch Wq and the branch W-i (see Fig. Q) correspond to this interval. 
Here one has to handle carefully with the interchange between branches Wo{x) and 
W-i{x) at the point x = — 1/e (this corresponds to the value a = 1). 
So, the nontrivial solution to Eq. (H) for the variable z is 
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(another choice of branches will be considered below). Thus, the solution to Eq. (g) 
we are interested in can be written in the form: 
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where we have introduced the function N(a) (see Fig. | 
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Fig. 3. The fiinrtinri N{a) (bold curve), and the functions No{a) (solid curve) and N—i{a) (dashed 
curve) (see Eq. (|l4[)). 

As follows from the definition (^) the function xe'-^ and the Lambert W function 
are mutually inverse. Here one has to distinguish precisely the branches Wq and 
W-i (see Fig. ^). From this figure it follows directly that for x > — 1 the function 
xe^ "corresponds" to the branch Wq, and for x < —1 the function a;e^ "corresponds" 
to the branch W-i. Let us introduce for convenience the notations of the inverse 



functions for these branches 



^Uwoy\x), x>-i, 

\(w.ir\x), x<-i. 



N{a) 



Now the function N{a) in Eq. ( |l4|) can be written in a more compact form 

-aWoUw^ir^i-a"^)], 0<a<l, 
-aW^i\{Woy\-~a-^) , a>l. 

It is worth noting that another choice of the branches of the function N{a) 

_ j iV_i(a), < a < 1, 
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leads to the trivial solution of Eq. (^ (in this case N{a) = 1). Therefore, the 
solution z ~ 1/N(a) does not satisfy Eq. (||) when a ^ I, and consequently we have 
to reject it. 

Let us briefly consider the basic properties of the function N{a) introduced in 
Eq. (|l^). In the physical range a > it is a nonnegative, monotone increasing 
function (see Fig. ^). With the help of the series (0)-(@) one can show that for the 
function N{a) the following expansions hold 
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4. Properties of the new analytic running coupling 

In this section we turn to the investigation of the properties of the one-loop new 
analytic running coupling in qcdI 

X^^(0^) = ^^, (21) 

For the convenience we omit the group factor and proceed with the expression 

a(Q2)^X^)(Q')-^- (22) 
z Inz 

Let us consider first of all the asymptotics of this function. In the ultraviolet (UV) 
limit the standard behavior of the invariant charge is reproduced: a{Q^ — s- oo) — s- 
1/lnz, z = Q^/A^, i.e., the asymptotic freedom of the theory is taken into ac- 
count. In the IR region there is an enhancement of the running coupling (|2^): 
a{Q^ 0+) — [elne]^^, e — Q^/A^. One should note here that such a behav- 
ior of the invariant charge is in agreement with the Schwinger-Dyson equationsj^ 



and as it was demonstrated recentlyjj'Q provides quark confinement without in- 
voking any additional assumptions. Thus the running couphng ( |2l| ) incorporates 
both the asymptotic freedom behavior and the IR enhancement in a single expres- 
sion. Obviously it is an essential advantage of our approach. It is easy to show 
that the function ( p^ ) has smooth behavior in the vicinity of the point = h?: 
aiQ"^ ^ A2) = 1 - S/2 + 0{S^), 5 = QVA^ - 1. It is worth noting also that for the 
running coupling ( p2|) the causal representation of the Kallen-Lehmann type holda3 
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In general, any expression for the running coupling makes sense only if the 
relevant definition of the parameter A is provided. Otherwise, a running coupling 
may be not a renorminvariant quantity at all. Thus, it is of a primary importance to 
present both the running coupling and the corresponding definition for its parameter 
A explicitly. 

So, let us represent the running coupling "a^jj in the renorminvariant form. In 
the general case for the invariant charge {Q'^ / , g) / (Air) = a{Q^) the normal- 
ization condition g'^{l,g) — g^ must be fulfiUed.Efl In our case this normalization 
condition acquires the form 



(mVA2) - 1 _Po 
ln(^7A2)/i2/A2 47r 
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Therefore, we have to resolve this equation with respect to the parameter A. The 
solution to such equations was considered in details in Sec. 3. With the help of the 
function N{a) (see Eq. dlj)) the solution to Eq. (|2^) can be presented in the form 



A^ = /i^ N 



(25) 



Thus, the renorminvariant expression for the one-loop new analytic running cou- 
pling (0) is the following 
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It is worth noting that when ^ oo (see Eq. (|T^) the right hand side of Eq. ( p5| ) 
tends to its standard formliJ corresponding to the perturbative running coupling 

<'(Q2) = 47r/(/3olnz) 
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Let us proceed to the construction of the [3 function corresponding to the running 
couphng (0). By definition,E^ 
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In the case under consideration, for the invariant charge ^al^^{Q'^) the (3 function 
takes the form 
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where z{a) is the solution of Eq. (||) with respect to the variable z (see Sec. 3). For 
the convenience we omit here the group factor, and use the expression "/San (a) = 
''/3an(a)/3o/(47r). With the help of the function N{a) defined in Eq. (||) the /3 
function for the NARC (UaI can be written in the form 
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Fig. 4. The solid curve presents the /3 funrtipn "/3an(a) (Eq. (|^)) corresponding to the one-loop 
new analytic running coupling '^'a^n (Eq. The standard perturbative result for this function 

/3s(a) = —a^ is also shown (dashed curve). 

Figure ^ presents this /3 function and the standard perturbative result /3s (a) = — 
corresponding to the one-loop perturbative running coupling a'^^Q'^) = 1/lnz. 
Thus, the /3 function corresponding to the one-loop new analytic running coupling 
"a^^jj reads as 

^a^-m-'i^^'^aiQ% %n(a) = J^f^. (31) 
Po z Inz Pq Po ln[A'(a)J 

The numerator and denominator on the right hand side of Eq. ( |30| ) have the 
opposite signs for all positive a. This follows directly from the expansion (^Sh (see 



below) and from Figs. || and ^. Thus, the (3 function ( |3C| ) possesses an important 
property, namely "/3an(a) < for all positive a. It is this property that provides the 
asymptotic freedom of the theory. By making use of the series (n8|)-(20) one can 



derive the following expansions of the [3 function (30). First of all 



%n(a -> 0+) 



O 



exp[ 

a 



(32) 



i.e., the standard perturbative limit is reproduced in the UV region. Secondly, 
%„(l + £) = -i-^e + 0(£'), e^O, 



(33) 



i.e., the function ( |30[ ) has a smooth behavior in the vicinity of the point a = 1. 
Next, in the IR region 
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Such a behavior of the (3 function provides the IR enhancement of the new ana- 
lytic running coupling. This enhancement leads ultimately to the quark-antiquark 
potential rising at large distances.ErQ 

Thus, as one could anticipate, the asymptotics of the obtained [3 function in- 
corporate both the asymptotic freedom and the IR enhancement of the running 
coupHng (0). 

Sometimes it is more convenient to present the (3 function in terms of the invari- 
ant charge g, a = g^/(47r). In this case the asymptotic expansions (|3|) and ( |3^ ) 
can be rewritten in the following way: "/3an(ff) — — /?o ^^/(Ifrzr^) when (this 
is nothing but the well-known standard perturbative resulttJ); and '^Pa.nig) — —g 
when g oo. It is worth noting here that a similar (i.e., linear in g) asymp- 
totic behavior (3{g) ~ —2g when^^ oo holds, when the so-called Mandelstam 
approximation is used (see papeiil3 and references therein). However the latter 
leads to the running coupling with more singular behavior in the IR region, namely 
a(Q2) ^ when ^ o+. 

As it was mentioned in the previous section, the new analytic running cou- 
pling (^ incorporates both the asymptotic freedom and the IR enhancement in 
a single expression. libviously this is an essential advantage of our approach. It 
was shown explicitlyO'El that invariant charge of the form (^ leads to the rising 
quark-antiquark (qq) potential without invoking any additional assumptions. The 
comparison of this qq potential with the so-called Cornell phenomenological poten- 
tial gives the following estimation for the Aqcd parameter: A = (0.60 ± 0.1) GeV. 
The fit has been performed in the physically meaning region0 0.1 fm < r < 1.0 fm 
(in this interval three active quarks should be taken into account, rif — 3). The anal- 
ogous comparison of the qq potential derived in Ref.B with the lattice dataEJ has 
been carried out recently. This gives a close value, A = (0.57 ± 0.1) GeV. Further- 
more, the estimation of the parameter Aqcd bas also been performed recently when 



deducing the value of the gluon condensate by making use of the NARC This 
gives the value A = (0.65 ±0.05) GeV. Therefore, we have the acceptable estimation 
of the parameter Aqcd in the framework of our approach: A — (0.60 ± 0.09) GeV. 

All this testifies that the new analytic running coupling (^ substantially involves 
the nonperturbative behavior of quantum chromodynamics. 

5. Conclusion 

The mathematical properties of the new analytic running coupling in QCD (|l|) 
are investigated. This running coupling naturally arises under analytization of the 
renormalization group equation^ One of the crucial points in our consideration is 
the explicit relation established between the NARC (Eq. (Q)) and its inverse function 
N(a) (Eq. ([T^)). The latter can be expressed in terms of the so-called Lambert W 
function. The properties of the function N{a) are examined in details. As known, 
when using the RG approach it is important to present the relevant quantities in a 
renorminvariant form. The function iV(a) enables one to perform this explicitly for 
the NARC (Q). Furthermore, the (i function corresponding to the running coupling 
"a^'j] is expressed in terms of the function N{a). The explicit form of the fimction 
"/^an in Eq. (^) provides an essential advantage when considering the gluon con- 
densate. The asymptotics of the /3 function in Eq. (^) are examined also. As one 
could anticipate, it incorporates both the asymptotic freedom and the IR enhance- 
ment of the invariant charge in a single expression. It was shown that there is 
a consistent estimation of the parameter Aqcd in the framework of our approach: 
A = (0.60 ± 0.09) GeV. Thus, the new analytic running coupling (|l|) substantially 
incorporates the nonperturbative behavior of the quantum chromodynamics. 

In further studies it would be interesting to investigate numerically higher-loop 
corrections to the NARC (|l|) and to compare the results with the explicit expressions 
obtained in the present paper. 
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